Enhanced Macroscopic Quantum Tunneling in Bi2Sr2CaCu208+5 
Intrinsic Josephson Junction Stacks 
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We have investigated macroscopic quantum tunneling in Bi2Sr2CaCu208+5 intrinsic Josephson 
junctions at millikelvin temperatures using microwave irradiation. Measurements show that the 
escape rate for uniformly switching stacks of N junctions is about A'^'^ times higher than that of 
a single junction having the same plasma frequency. We argue that this gigantic enhancement of 
macroscopic quantum tunneling rate in stacks is boosted by current fluctuations which occur in the 
series array of junctions loaded by the impedance of the environment. 

PACS numbers: 74.72.Hs, 73.23.-b, 73.40.Gk, 85.25.Cp 



After 20 years since its discovery, macroscopic quan- 
tum tunneling (MQT) in Josephson junctions remains 
a fascinating phenomenon which attracts interest of a 
broad physics community. MQT was first observed in 
Nb Josephson junctions at very low temperatures ,1] and 
has been used to study energy level quantization by mi- 
crowave absorption More recently, so-called phase 
qubits based on MQT in current-biased Josephson junc- 
tions have been reported as very promising hardware for 
quantum information processing . 

Due to their rf-wave order parameter symmetry 
cuprate high-T^ superconductors were initially regarded 
unsuitable for MQT experiments. However, it has been 
argued that although there are nodes in the d-wave order 
parameter, MQT should not be suppressed completely in 
high-Tc Josephson junctions Meanwhile, recent mea- 
surements on YBa2Cu307_5 grain boundary junctions 
show the solid evidence of MQT 's'l . 

On the other hand, intrinsic Josephson junctions (IJJs) 
in layered high-Tc superconductors |^ ^3 are rather 
attractive candidates for MQT experiments. IJJs in 
Bi2Sr2CaCu208+5 are formed by pairs of Cu02 dou- 
ble planes, separated by a Bi2 03 insulating layer. They 
have a much higher Josephson coupling energy than grain 
boundary junctions and a better homogeneity, as they are 
located inside a more or less perfect single crystal. More- 
over, intrinsic junctions exhibit current transport along 
the c-axis direction, i.e. perpendicular to the copper ox- 
ide layers, where the nodes of the dx2-y2 order parameter 
should not affect MQT at all. With the recent invention 
of a double-sided fabrication technique i O'^s can avoid 
heating the junctions by contact resistance. This makes 
MQT in IJJs practically attainable. Recently, MQT has 
been observed on a single junction in Bi2Sr2CaCu2 08-i-(5 
IJJs [Ulli. The temperature T* of the crossover be- 
tween thermal and quantum escape was reported to be 
rather high. 

In this Letter, we present an experimental study of 
MQT in IJJs using microwave spectroscopy. We demon- 
strate that the unique uniform array structure of intrinsic 



Josephson junction stacks causes an enormous enhance- 
ment of the tunneling rate. We argue that this enhance- 
ment can be caused by current fluctuations in the stack. 

The samples were fabricated using the standard 
double-sided ion beam etching technique |^. The stack 
height ranges from approximately 7.5 nm to 15 nm, i.e. 
the number of junctions N in our IJJ series arrays is be- 
tween 50 and 100. The junction area varied for differ- 
ent stacks from 1x2 /im^ to 2x3 fim^. Critical current 
densities jc were between 0.7 - 2.8 kA/cm^. Sample pa- 
rameters are summarized in Table 1. Current-voltage 
(I-V) characteristics were measured with a standard four- 
terminal configuration using current biasing. Experi- 
ments were performed in a "^He cryostat with a minimum 
temperature of 300 mK and in a ■^He/^'He dilution refrig- 
erator with a base temperature of 10 mK. 

We investigated two different types of samples. Their 
typical I-V characteristics are presented in Fig. 1. In the 
first case, at least one junction had a critical current sig- 
nificantly lower than those of the rest of the stack. This 
was achieved either by exposing the unprotected side of 
the crystal to ambient atmosphere, resulting in a reduced 
critical current density, or by over-etching resulting in a 
reduced cross-sectional area. When increasing the cur- 
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FIG. 1: I-V characteristics: (a) One-by-one switching of three 
junctions in a stack, ttSJ3; (b) uniform switching of the whole 
stack, (jUSl. The insets show the full-range I-V curves for 
both samples. The sample geometry is sketched in the upper 
right corner. 
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FIG. 2; Temperature dependence of the standard deviation 
of the switching current from the zero-voltage state. Solid 
circles: sample tiUSl; open circles: sample USJl. The straight 
line is a linear fit to all data points above 700 mK. The inset 
shows switching current distributions for [tUSl. For better 
comparison, the IJUSI data were divided by a factor of 1.09. 



rent, we observed voltage jumps from one branch to the 
next in steps of approximately 25 mV (Fig. la). Due 
to the large difference of the critical current, it is always 
the same junction that switches to its resistive state first. 
We call this behavior "single-junction switching" and de- 
note this type of samples as the jJSJ series. In the second 
case, the stack is so homogeneous that upon increasing 
the current, we never observed any stable states between 
zero voltage and 1.14 V, where all junctions are in the 
resistive state (Fig. lb) . We call that behavior " uniform- 
stack switching" and denote this type of samples as the 
jJUS series. Nevertheless, we find that many branches 
corresponding to different numbers of resistive junctions 
are still there. They can be easily traced up from the 
return curve of the resistive state |eIIioIIii|. 

Distributions of the switching current, Is, at which the 
junctions escape from the zero voltage state [3| were 
measured using a high-resolution ramp-time based setup 
|l5l |. The bias current was ramped up with a rate of 
200 mA/s. After detecting a switching event by a volt- 
age threshold of 20 /iV, the current was switched to zero 
within less than 10 /its in order to avoid heating effects. 



TABLE I: Sample parameters measured. 
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Monitoring the voltage on a fast oscilloscope for jJUS- 
type samples showed that the current decrease was fast 
enough so that any stack always switched to its first re- 
sistive branch, i.e. to a state with one resistive junction. 
Based on this fact we emphasize that regarding dissipa- 
tion both types of samples were measured under exactly 
the same experimental conditions. Switching current dis- 
tributions were measured at temperatures between 20 
mK and 12 K using a repetition rate of 600 Hz. The 
switching current statistics was determined from 20,000 
to 60,000 switching events. The switching probability 
P(I) shown in the inset of Fig. 2 is defined as the number 
of switching events per /iA normalized to the total num- 
ber of events. The standard deviation a of the switching 
current was determined from the width of the P(I) curve. 
Fig. 2 shows cr as a function of temperature for samples 
(jSJl and ttUSl. The saturation of o'(T) at low tempera- 
tures corresponds to a crossover from thermal activation 
to MQT. For sample jjSJl the saturation is not complete 
in Fig. 2. This is due to the fact that this experiment was 
done in a ^He cryostat where the lowest temperature was 
300 mK. We verified complete saturation in our experi- 
ments in the dilution refrigerator. Nevertheless, one may 
still see in Fig. 2 that the crossover temperature T* of 
the single-junction sample ftSJl is about 300 mK, while 
the uniform- switching sample jJUSl shows T* of about 
700 mK. This discrepancy can not be fully accounted by 
the difference in the critical currents. 

At temperatures below T*, we measured the switching 
current distributions under microwave radiation in the 
frequency range between 10 and 40 GHz. Such measure- 
ments allow to determine the plasma frequency ujp and 
the absolute value of the fluctuation-free critical current 
Ic directly Q . Here, microwave spectroscopy serves as a 
tool to determine the energy level separation in the quan- 
tum regime. It should be noticed that quantum transi- 
tions between levels cannot be easily distinguished from 




FIG. 3: (a) Density plot of switching current vs. microwave 
power for two-photon absorption in sample JUS3 at 30 mK. 
The switching probability P(I) is color coded according to the 
bar on the right-hand side, (b) The corresponding enhance- 
ment of the escape rate at different microwave powers. The 
inset shows the fit of the 6 dBm enhancement curve with a 
Lorentzian. 
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FIG. 4: Applied microwave frequency versus normalized res- 
onance current for sample jiSJ2 (open circles) and ttUS3 (solid 
circles) measured at 20 mK. The data are fitted to Eq.(l) 
using n=2 (upper dotted line and upper solid line) and n=3 
(lower lines) . The vertical lines mark the maximum switching 
current of the samples. 



classical plasma resonance peaks using just spectroscopy 
data [l^. However, below the crossover temperature T* 
quantum fluctuations dominate thermal ones and quan- 
tum mechanics is appropriate for describing the system. 
Since the zero-bias plasma frequency = Wp/27r for the 
samples was well above 100 GHz, we used multi-photon 
|l7l | rather than single-photon absorption in our measure- 
ments. 

As the microwave power is increased, the P(I) distri- 
bution becomes double-peaked. The microwave-induced 
peak at lower currents corresponds to a plasma resonance 
in the junction. In the quantum picture, this peak is 
interpreted as tunneling from highly-populated energy 
level(s). Figure 3(a) shows a density plot of the switch- 
ing current distribution of a uniform-switching stack ver- 
sus microwave power at 38.2 GHz. Figure 3(b) shows 
the corresponding enhancement of the escape rate. The 
double-peaked P(I) distribution develops at a microwave 
power of about -1 dBm referred to the top of the cryo- 
stat. The two-photon resonance current peak appears 
at about 57.0 /iA. It should be noted that the width of 
the resonance peak is smaller than the distribution width 
at zero microwave power, which indicates that the stan- 
dard deviation of P(I) measured without microwaves is 
not limited by current noise in our setup. The resonance 
current is defined as the position of the resonant peak 
when both peaks have equal amplitudes. The results at 
different microwave frequencies are summarized in Fig. 4. 
The data are fitted to 
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(1) 



where ojp is the plasma frequency of the junction at the 
normalized bias current 7, is the plasma frequency at 
zero bias, and n is the number of photons taking part 
in the absorption process. 7 = Ib/Ic is given in nor- 



malized units, where lb is the bias current and Ic is the 
fitted fluctuation-free critical current. The data in Fig. 4 
show the best fit to Eq. (1) by assuming two- and three- 
photon absorption. At lower frequencies and high pow- 
ers we observed multi-photon peaks up to n = 6 (not 
shown). From the fits we obtain = ^p/27r = 150 
GHz (sample jJUSS), and = 180 GHz (sample tJSJ2). 
The fitted /p of the other samples are shown in Table 1 . 
The obtained plasma frequencies provide an estimate for 
the junction capacitance per unit area C=jc/(27r$o/p)i 
where $0 is the magnetic flux quantum. Our measure- 
ments yield G=3.9 /xF/cm^ and £^=5.3, which conforms 
well with the value of £^=5 obtained in earlier work p^ . 

From the results shown in Fig. 4 and Table 1 , one may 
notice that the normalized switching current 75. = Is/ Ic 
of the uniformly switching stacks is significantly lower 
than that of the single junction switching samples. The 
typical value of 7^ for single-junction samples is about 
99% of the fiuctuation-free critical current Ic, whereas 
for uniformly-switching stacks it is only 96%. The switch- 
ing range of the single-junction samples is fairly similar 
to that of Nb junctions with comparable critical current 
densities, where 7s is also close to 99% On the 

contrary, the normalized switching bias current 7s of uni- 
formly switching stack samples is about 3% smaller. This 
discrepancy can explain the difference in T* for [tSJ and 
ttUS samples. The plasma frequency lo^ at the switch- 
ing current is given by Eq. (1) with 7 = 75. For a sin- 
gle junction, the temperature of the crossover from ther- 
mally activated escape to MQT is expected 19] to be 
T* K, hujp/2TrkB- For uniformly-switching samples, the 
premature switching at lower 7s leads to the high uj^ and 
thus to an enhanced T*. 

In the quantum regime, a single Josephson junction 
escapes through the energy barrier with an escape rate F 
determined by 



^ ojp /864C/o7ry/' f 36 Uo 
i = I — ; I exp 



27r \ hujp 



where Up is given by Eq. (1), Uq = Ej^^{l - 7)^/^ 
is the height of the potential barrier at 7 < 1, and 
Ej = ^o-^c/27r is the Josephson coupling energy. In Fig. 5 
we plot Eq. (2) using the values of ojp and Ic in Table 1 for 
samples jJUSl, ttUS4 and ttSJ3, respectively (solid fines). 
The interconnected dots were determined from the mea- 
sured P(I) histograms. For single junction sample ttSJ3, 
one can see that experiment agrees well with the theory. 
However, there is a huge difference between theory and 
experiment for both uniformly switching samples. 

We used the mean switching bias 7s of each sample to 
mark the intersection of 7s with the theoretical curves in 
Fig. 5 by large open dots. According to the theory, these 
dots should correspond to the most probable escape cur- 
rent. One can see that the intersections of 7s with the 
experimental curves show actual escape rates higher by 
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FIG. 5: Escape rate versus bias current for uniformly switch- 
ing stack samples (tUSl and ttUS4 and single-junction switch- 
ing sample (jSJ3. The solid lines are calculated by the single- 
junction quantum escape theory (2). The points intercon- 
nected by lines are extracted from the measured P(I) distri- 
butions. The vertical-arrows indicate the enhancement of the 
escape rate (see text). 

factor of roughly 2000 (tJUSl) and 10000 (ttUS4). This dif- 
ference is marked by the vertical arrows in Fig. 5. Sam- 
ples listed in Table 1 have different sizes and j^s. We 
have measured them in different setups and at different 
temperatures. Experiments for all single junction sam- 
ples well agree with theory (2). On the other hand, for 
all uniformly switching samples we always found a huge 
difference between the measured data and single-junction 
model given by Eq. (2). We therefore conclude that the 
single junction tunneling theory is not applicable to the 
uniform switching stacks. 

The value of N in jJUS samples can be determined 
by counting the number of resistive branches in the I-V 
characteristics. Assuming that the single junction escape 
rate is enhanced by a factor of , one can renormalize 
the peak escape rate of all jjUS samples. In Fig. 5 we 
have chosen to show samples jlUSl with = 46 and 
jiUS4 with TV = 100. We find that the corrected val- 
ues shown by large solid dots conform pretty well to the 
experimental data. Moreover, our data for samples ttUS2 
and jlUS3 (see Table I) also match this N'^ correction. 

If there would be no interaction between TV identical 
junctions in a current-biased series array, the escape rate 
F (the probability that at least one junction switches at 
a given bias current) should be merely enhanced by a 
factor N with respect to a single junction. If we assume 
that switching of any junction in the array is triggered, 
in addition to its own fluctuations, by its nearest neigh- 
bors in the array (e.g., via charge coupling through the 
shared Cu02 double-planes), then we should get an en- 
hancement by factor of about 3A'^. However, experiments 
imply that the enhancement of the escape rate in stacks 
is proportional to 7V^. This is possible when there is 
an interaction between any pair out of the N junctions. 



Such interaction occurs when the stack of junctions con- 
nected in series is loaded in parallel by a relatively low 
impedance. Fluctuations of the phase difference of a sin- 
gle junction change its parametric Josephson inductance 
and thus the total inductance of the array. The external 
bias current is split-up between the array and the exter- 
nal impedance, which can be regarded as the impedance 
of the environment at the plasma frequency. The bias 
current flowing through the array thus changes under 
fluctuations in any junction. A speciflc analysis of this 
model goes beyond the scope of this experimental paper. 

In conclusion, we have found a drastically enhanced 
escape rate of macroscopic quantum tunneling in uni- 
formly switching Bi2Sr2CaCu208+5 intrinsic Josephson 
junction stacks. This enhancement adds a factor of ap- 
proximately N'^ to the quantum escape rate of a single 
Josephson junction. This can be caused by large quan- 
tum fluctuations due to interactions among the N junc- 
tions and results in a significant increase of the crossover 
temperature T* between the thermal activation regime 
and quantum tunneling. 
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